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^N) ■ Abstract. Equilibria of polygonal linkage with respect to Coulomb potential of point 

t, ' charges placed at the vertices of linkage are considered. It is proved that any convex 

O ' configuration of a quadrilateral linkage is the point of global minimum of Coulomb po- 

^(^ , tential for appropriate values of charges of vertices. Similar problems are treated for the 

^^ 1 equilateral pentagonal linkage. Some corollaries and applications in the spirit of control 

f"C^ ' theory are also presented. 
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1. Introduction 

Preliminary remarks. Motivated by the famous Maxwell conjecture on equilibria of point 
charges ^ (cf. also ^) we deal with the the Coulomb potential of a system of point charges 

C^ I placed at the vertices of a (flexible) planar polygonal linkage. We consider Coulomb potential 

of a vertex-charged linkage as a meromorphic function on the planar moduli space of the 
linkage and investigate its critical points. 

The scenario we have in mind is suggested by some recent research concerned with the 

^ I control of nanosystems [?]■ As an abstract analog of a real physical situation we suggest 

^^ ' the following setting which can be described as controlling the shape of linkage by the values 

of charges at its vertices. The basic implicit assumption and motivation is that a vertex- 
charged linkage subject only to Coulomb interaction of its charged vertices should take the 
shape with the minimal Coulomb potential. 
'nI . This setting suggests several aspects and problems. In the present paper, we concentrate 

on the following scenario. Civen a planar configuration of linkage we wish to find the 
vertex charges such that the global minimum of the arising Coulomb potential is achieved 
at the given configuration. Such a collection of charges will be said to stabilize the given 
configuration. If any configuration of linkage has a stabilizing system of charges we will say 
that this linkage admits a complete Coulomb control. 
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H ■ Assuming that a stabilizing collection of charges exists what does the set of all stabilizing 

. .' charges look like? How many minima and other critical points has a system of stabilizing 

charges? 

An interesting special case of the previous problem arises if one asks if any convex con- 
figuration can be stabilized by a system of charges of the same sign. This problem will be 
called Coulomb control of convex configurations. 

It should be noted that this research arose as a natural continuation of our previous joint 
results on Morse functions on moduli spaces of polygonal linkages [4] , [5] , [6] . The present 
paper was completed during a " Research in Pairs" session in CIRM (Luminy) in January of 
2013. The authors acknowledge the hospitality and excellent working conditions at CIRM. 

Definitions and results. A polygonal linkage L is defined by a collection of positive num- 
bers I = (/i, ...,ln), called sidelengths, which we express by writing L — L(l). 



Key words and phrases. Polygonal linkage, planar configuration, point charge, Coulomb potential, critical 
point, equilibrium. 



2 G.KHIMSHIASHVILI*, G.PANINAt, D.SIERSMA* 

Physically, a polygonal linkage is a collection of rigid bars of lengths k joined in a cycle 
by revolving joints. It is a flexible mechanism which can admit different shapes, with or 
without intersections. 

By M{L) we denote the moduli space of planar configurations, that is, the space of all 
polygons with the prescribed edge lengths factorized by isometrics of R^ : 

M{L) =: {(pi, ...,Pn)\pi e M^, \piPi+l\ = k, \pnPl\ = ln}/Iso{M.^). 

This is not exactly the moduli space M-{L) treated in 4 and 2], where the space of 
polygons is factorized by orientation preserving isometrics. However, there is a two-fold 
covering M{L) -^ M{L). 

By M^ {L) we denote the set of all convex configurations. We allow here non-strictly 
convex polygons, that is, those having (at least) one angle equal to tt. The latters obviously 
form the boundary dM'^'{L). The set of all strictly convex configurations (all angles are 
less than tt) is the interior IntM'^ . It is known (see [9]) that M'^'{L) is homeomorphic to 
a ball. In this paper we only deal with n — 4, 5. For a 4-bar polygonal linkage, M{L) is a 
(topological) circle, whereas M'~'{L) is homeomorphic to a segment. For a 5-bar polygonal 
linkage, M{L) is (generically) a surface, whereas M'~^ [L) is the disk D^ . 

Putting a collection of charges qi at the vertices pi of a configuration and considering the 
Coulomb potential of these charges we get a function defined on M{L). We will refer to this 
setting by speaking of a vertex- charged linkage with the system of charges qt. 

Recall that the Coulomb potential £■ of a system of point charges g^ £ K placed at the 
points Pi of Euclidean plane R^ is defined as 

^-Ef^' (1) 

where Xij = \pi — pj\ is the distance between ith and jth charges. 

Since we are only interested in critical points of Coulomb potential, addition of a constant 
makes no difference. By the very definition of polygonal linkage, the distances corresponding 
to two consequent vertices in formula ^ remain the same for all configurations of linkage. 
Hence their sum is constant for any fixed collection of charges and for our purposes it is 
sufficient to work with the effective Coulomb potential E of configuration defined as 

^ = E!^' (2) 

where Xij is the length of diagonal between (non-neighboring) ith and jth vertex of the 
configuration. We say that a collection of charges stabilizes the configuration P \i E attains 
at P its minimal value. In this case we say that P is the minimum point of E. 

We explicate now the setting and notation. For n = 4, we put one positive charge t at 
the first vertex. The rest three vertex-charges are -1-1. 

For n = 5, we put two positive charges s and t at any two non-neighboring vertices 
and say that s and t are controlling charges. The rest three charges, called non- controlled 
charges^ are again equal to -1-1. 

Remark 1. // all the pairs of consecutive sidelengths are different, E is a smooth function 
without poles. If not, in our setting we have only the poles with positive residues, which do 
not affect our study of minima points of E. 
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Our main results are: 

(1) For any 4-bar linkage, and any positive charge t, there exists a unique convex con- 
figuration which is a critical point of E. This is the global minimum of E, and it 
depends continuously on t. 

(2) For any 4-bar linkage, and each convex configuration P, there exists a unique charge 
t which (together with the non-controlled charges) stabilize P. In this case, t is 
positive, and P is the global minimum point. 

(3) We have a complete Coulomb control for convex quadrilaterals. This means a two- 
step navigating algorithm bringing any 4-bar configuration to an in advance pre- 
scribed convex position ruling by a positive charge t. 

(4) For any convex equilateral pentagon P, there exists a unique pair of positive charges 
(s, t) for which P is a critical point of E. However, it is unclear whether P is the 
global minimum point. 

2. Coulomb control problem of convex quadrilaterals 

We begin with considering Coulomb control for convex configurations of a non-degenerate 
4-bar linkage L with one positive controlling charge t > at the first vertex and three equal 
charges +1 at the other three vertices. For a convex configuration P of L, we denote by 
X and y the lengths of its two diagonals, and by E its (effective) Coulomb potential with 
controlling charge t: 



^ t 1 

E=- + -. 
X y 



(3) 



Lemma 1. For a given convex quadrilateral P G M{L), there exists a unique t such that P 
is a critical point of E on M{L). In this case, t is positive. 

Proof. In a neighborhood of P in M{L), we have a relation of the form y = y{x). Then 
the condition that P is a critical point of E is 

y'^/x^ == -ty'{x), 

which defines t uniquely. For a convex configuration P, each flex of P which increases one of 
diagonals, shortens the other one. This means that y' {x) is negative. Hence t is positive. D 

Proposition 1. For a 4-bar linkage and a positive charge t, 

(1) E has a unique minimum in the interior of M'~^{L) (that is, among strictly convex 
configurations), which depends on t continuously. 

(2) E has no critical points among non-convex non-intersecting 4-gons. 

(3) There is at least one critical point of E, which is a self-intersecting ^-Qon. 

Proof. 

(1) First, we show that all the points (x^,y^) lie on an elliptic curve. 

Indeed, let 
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Put 

F{x, y) = det(M(x, y)) = -2x^y^ - 2x^y'^ + 2x^{d^ - a^){b'^ - c"^) + 2y^{b^ - a^){d^ - c"^) 

-2(ac - bd){ac + bd){a^ -h^ ^c^ - d^). 



G.KHIMSHIASHVILI*, G.PANINAt, D.SIERSMA 



it 



which is the Cayley-Menger determinant. The equation F{x, y) = defines a real algebraic 
curve of degree 6, which depends only on x'^ and y^. In fact it is a cubic g{w,z) with 
w = x'^, z = y^ . The configuration space M{L) maps bijectively to the compact convex 
component {oval) of the curve (inside a box bounded by squared maximal length of diago- 
nals), see Fig. [1] This component has has not more than two intersection points with any 
straight line. 
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Figure 1. Diagonal relation in x,y (left) and in x'^,y'^ (right) 




Figure 2. The configuration space M{L) is divided into four parts 



Consequently, the oval contains exactly 2 points where §f- = and two other points 
where ^ = 0. These points divide the oval into 4 segments, each determined by the signs 
of partial derivatives (see Fig. [5]). 
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(1) Convex 4-gons (that is, M'-^) correspond to the " ++" case. A vertex of configuration 
is aligned if the angle adjacent to it is tt. For a non-degenerate 4-bar hnkage, there 
exist exactly 2 vertices which can be aligned. The corresponding configurations are 
the boundary points of AI^' . 

(2) Non-convex and non-self-intersecting 4-gons correspond to "H — " and " — h" parts. 

(3) Self-intersecting 4-gons correspond to the " " part. 

We prove now that, for a convex P, we have y"^ < 0. 

Indeed, locally m; is a function in z. We have w = 10(2), therefore y'^ — w{x'^), which 
implies 

y = w^^Hx^) 
I -1/2 / 

We know from convexity of the elliptic curve that u'"^ < 0. 

y^^ U;2 u;3/2 + y,l/2 '^ "■ 



Assume that P is a convex configuration which is critical for E. 

t I . ,. 
E — — \ — miplies 
x y 

11' 

and 





F' -* y'^ 

" x^ y2 ' 




_ 2i ^ 2[y'{x)]^ y"{x) 

X^ y3 y2 



which is positive. 

Therefore, each critical point in the convex part is a local minimum. 
So the statement (1) of the proposition is now straightforward: if there are two local 
minima, there should be a local maximum in between, which is impossible. 

(2) For a non-intersecting 4-gon there is a flex strictly increasing both diagonals. 

(3) There should be a maximum point of E in M{L). By what is proven above, the 
domain of self-intersecting 4-gons is the only possibility for it. D 

Theorem 1. (1) We have a continuous hijection 

^ : (0,00) -^ IntM'^iL) 

which maps a charge t to the minimum point of E. 
(2) The mapping Lp extends to a hijection Lp : [0, 00] — > M (V). 

Proof. Point (1) follows from Proposition (TJ (1) and Lemma[T] Point (2) follows from (1) 
and the above remark. D 



Lemma 2. Assume that L is a J^-har linkage, and the charge t is zero. Then E has exactly 
two critical points on M{L): one minimum and one maximum. 

Proof. The problem is reduced to finding the critical values of the diagonal x. D 

The lemma means that if we start at any configuration of any 4-gon L and put zero charge 
at the ruling vertex, the gradient flow will bring us to the global minimum. 
However, for other values of i, there might be extra local minima: 
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Example 1. For a = 6, b = 6.5, c = 6.2, d = 5.8 and t — 2, numerical computations show 
that there are four critical points of E in total (see Fig. \3^. The lengths of the diagonals 
and the values of E (up to two decimals) are: 

(1) X = 0.50, y = 3.24, E = 2.61 (local minimum) 



X 

(2) X 

(3) X 

(4) X 



' 4.11, y = 0.30, E = 6.90 (maximum) 

= 1.24, y = 0.58, E = 4.24 (local maximum) 

- 9.59, y = 7.60, E = 0.36 (minimum) 




x=9.59; y=7.60 
absolute minimum 
E=0.36 



x=1.24;y=0.58 
maximum 
E=4.25 




B C 




D 



x=4.11;y=0.30 
maximum 
E=6.9 



B 

x=0.50; y=3.26 

local minimum 
E=2.61 



Figure 3. Critical configurations from Example [T] 

The above leads us to the following navigating algorithm. 

Algorithm 1. Assume we are given a 4-bar linkage L, its unknown starting configuration 
P and a target convex configuration P' . The following algorithm tells us how to reach P' by 
altering the ruling charge t. 

(1) Put the ruling charge t = 0. 

(2) Compute the stabilizing charge t' of the configuration P' . 

(3) Put the ruling charge t = t' . 



Remark 2. If we allow non-positive charge and aim at non- convex configurations, the 
situation becomes more complicated. Example [7] shows that the number of (local) minima is 
greater than 1. 

Remark 3. If we put the same charged 4-bar linkage in 3D, we can skip step 1 in the above 
algorithm. Indeed, all critical points of E are obviously planar configurations. Besides, 
unlike the 2D case, a local minimum cannot be self -intersecting: the unfolding flex increases 
both of the diagonals. Therefore, in 3D the potential E has exactly one local minimum which 
is the global minimum. 

Remark 4. All above lemmata and Theorem 1 remain valid if we replace the Coulomb 
potential E either by 



E- = E 



[Xij ) 



a > 0, 



(4) 



or by the limit version of (Rl) 



i?'"=E9'9jln(^ 



(5) 
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3. Coulomb control of convex equilateral pentagons 

Let L be a equilateral 5-bar linkage. In our setting we put two positive charges s and 
t at a pair of non-neighbouring vertices, say, at the third and fifth vertex of configuration. 
The charges of other vertices are constant and equal to +1. Then the (effective) Coulomb 
potential E : M{L) — > M is defined by the formula 



E = 



1 



1 



t 



st 



xi4 X2A a;i3 2:25 



X3b 



(6) 



where Xij are the lengths of diagonals of the configuration V . 

We wish to understand whether two non-adjacent charges can provide a full control (in 
the same sense as for quadrilateral linkages) on convex configurations of L. 




Figure 4. 

Lemma 3. Let P G dM'-^ he a critical point of E for some charges {s,t). Then {s,t) never 
belongs to (0, 00) x (0,oo). 

Proof. Assume the contrary: the charges s and t are both from (0, 00). A configuration 
is a boundary point of AI^{L) whenever one of the vertices is aligned. Then by pulling the 
aligned vertex orthogonally to the adjacent edges (see Fig. [S]) we get an infinitesimal flex 
which yields a flrst-order decrease of E. O 




Figure 5. 



Lemma 4. For s < 0, i > 0, a convex pentagon is never a critical point of E. 
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Proof. Assume the contrary. Fix X13 and compress X25, see Figure El This yields an 
infinitesimal flex such that E' < 0. D 




Figure 6. 

Theorem 2. For each strictly convex configuration P of an equilateral pentagonal linkage, 
there exists exactly one {s,t) £ (0,oo) x (0, 00) such that P is a critical point of E for the 
charges {s,t). 

Proof. We start by reminding that 

E = f /xi4 + 1/2:24 + t/xi3 + s/x25 + St/X35. 

Take the diagonals X35 and X13 as local coordinates in a neighborhood of P. 
The polygon P is a critical point of E means that dE vanishes: 

st 



'dE/dx35 = ai/a;?4 + /3i/x24 + sji/xl^ + ^- = 0, 



and 



f-ss 



-dE/dxi3 = a2/xi4 + h/x'^i + t/xfg + ^ = 0, 



'-25 



where 



ai = 8x14/8x35, Pi = 8x24/8x35, 71 = 8x25/8x35, and 



a2 = 9x14/9x13, /3i = 8x24/8x13, 71 = 8x25/8x13 
We get a system in two variables s and t which reduces to the following quadratic equation 
in s: 

A + Bs + Cs^ = 
with 

^2 ^ ^2 ' 
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'''25 -^35 •''14 •'•24 



and C 



xhl2 

2 2 ■ 
■''35^25 



It is sufficient to prove that AC is negative. Then the equation has exactly one real 
positive solution s. By Lemma SI t is also positive. 

Indeed, AC is negative (and the proof is completed) because of the following sign analysis. 

(1) Consider the flex of P which fixes Xi^ and stretches x^^. 

Xii compresses, therefore ai = dxu/dx^c, < 0. 
X2i compresses, therefore /3i = 8x24/0x3^ < 0. 
2:25 compresses, therefore 71 = 8x25/8x35 < 0. 

(2) Similarly, fixing 0:35 and stretching 113, we obtain that a2 < 0, /32 < 0, and 72 < 0. 
Let us give a more detailed explanation why the above inequalities hold. As an example, 

let us prove that if the diagonal X13 is fixed, and the diagonal X35 gets longer, then a;24 
compresses (see Fig. [7]). 

We can assume that the three points 1,2,3 are fixed and that the points 4 and 5 are 
moving with infinitesimal velocities V4 and V5. The vector v^ is orthogonal to the edge 43. 
By convexity, the angle between V4 and the diagonal 24 is smaller than 7r/2, which implies 
the claim. D 




Figure 7. 

Remark 5. If we try to rule by putting charges at adjacent vertices (say, at the vertices 
1 and 2), the situation gets worse: we reach not all of the convex polygons. For instance, 
we will never he able to have simultaneously the vertices 5 and 3 aligned. By continuity, an 
entire neighborhood of such a polygon becomes unreachable. 

4. Concluding remarks 

Our results suggest several natural problems and conjectures concerned with the Coulomb 
control scenario for n-bar linkages with arbitrary n. In particular, there is good evidence 
that, for equilateral linkages, the complete Coulomb control of convex configurations is valid 
if we are permitted to choose controlling charges at ALL vertices. 
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Assuming that this conjecture is true, what is the minimal number of controlling charges 
sufficient to provide complete Coulomb control? For dimension reasons, one may await that 
it should be sufficient to have n — 3 controlling charges. Is this amount of controlling charges 
indeed sufficient for any n-bar linkage? 

Several interesting problems in the spirit of the famous Maxwell conjecture [8] are con- 
cerned with estimating the amount and possible types of critical points of Coulomb potential 
in the moduli space of linkage. In particular, what is the exact upper bound for the number 
of critical points of Coulomb potential in the planar moduli space of generic n-bar linkage? 
We conjecture that, for a generic 4-bar linkage, the Coulomb potential can have not more 
than four critical points and their types are the same as in Example 1. 

Analogous problems are also meaningful and interesting for all potentials of the form 
O, ([5|). These and further aspects of the Coulomb control scenario will be addressed in 
forthcoming publications of the authors. 
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